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III. " On Jacobi's Figure of Equilibrium for a Rotating Mass of 
Fluid." By G. H. Darwin, M.A., LL.D., F.R.S., Fellow of 
Trinity College and Plumian Professor in the University of 
Cambridge. Received October 12, 1886. 

I am not aware that any numerical values have ever been de- 
termined for the axes of the ellipsoids, which are figures of equi- 
librium of a rotating mass of fluid. * 

In the following paper the problem is treated from the point of 
view necessary for reducing the formulae to a condition for computa- 
tion, and a table of numerical results is added. 

Let a, b } c be the semi-axes of a homogeneous ellipsoid of unit 
density ; let the origin be at the centre and the axes of x, y, z be in 
the directions a, b, c. 

Then if we put — 

A*~a? + u, J3 3 =& 2 + ^, C 2 =c 3 +^, and 

r°° du 

*=) lBO <!> 

it is knownf that the potential of the ellipsoid at an internal point 
x, y, z is given by — 



TT 7 T T , %? d^ , y%d& z* d¥ 

Y=urabo\ * + — T- + 1-37-+— -7- 
L a da do c do m 

Now let us introduce a new notation, and let 



• • • * • \^ ) 



* The following list of papers bearing on this subject is principally taken from a 
report to the British Association, 1882, by W. M. Hicks : — 

Jacobi, ' Acad, des Sciences/ 1834 ; Liouville, ' Jonrn. !Ecole Polytech./ vol. xiv, 
p. 289 ; Ivory, ' Phil. Trans.,' 1838, Pt. I, p. 57 ; Pontecoulant, ' Syst. du Monde,' 
vol. ii. The preceding are proofs of the theorem, and in more detail we have : — 
0. O. Meyer, 'Crelle,' vol. xxiv, p. 44; Liouville, ' Liouville's Journ.,' vol. xvi, 
p. 241 ; a remarkable paper by Diriehlet and Dedekind, * Borchardt's Journ.,' 
vol. lviii, pp. 181 and 217 ; Kiemann, ' Abh. K. Gres. Wiss. Grottingen,' vol. ix, 1860, 
p. 3 j Brioschi, ' Borchardt's Jonrn./ vol. lix, p. 63 ; Padova, ' Ann. della Sc. 
Norm. Pisa,' 1868-9 (being Diriehlet and Biemann's work with additions) ; 
G-reenhill, * Proc. Camb. Phil. Soc.,' vol. iii, p. 233 and vol. iv, p. 4 ; Lipschitz, 
< Borch. Jonrn.,' vol. Ixxviii, p. 245 ; Hagen, ' Schlomilch Zeitsch. Math.,' vol. xxiv, 
p. 104 ; Betti, < Ann. di Matem.' vol. x, p. 173 (1881) ; Thomson and Tait's ' Nat. 
Phil.' (1883), Part II, §778 ; a very important paper by Poincare, * Acta Mathem.,' 
7, 3 and 4 (1885). 

f Thomson and Tait's 'Nat. Phil.' (1883) §494, I. The form in which the 
formula is here given is slightly different from that in (8), (11), (15) of §§ 494, k, I. 
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c=acos7, sma = 




a 



.»— 6* 



a 3 ~~c 3 ' 



and b= a cos ft, 



(3) 



so that sin j3 = sin « sin 7, and & = a s/ (1 — sin 3 * sm 3 <y .) 



»3 — _„ /»3 om*(' 



Also let jL 2 =w + a 2 =—T-7rr= & 3 -^-37;' 



1 



whence 



and 



and 



3 =^ + c 3 =^$-cos^ 

, 2a 3 sin 2 7 „ , rt 

d^= — — r-57— cos d9, 

surO 



Y 



(4) 



/•« 



c?%=2a%in 3 7 







cos0 , 09.9 f Y cos7 . 
-t-^-: a0=2a 3 sin 3 7 —-J-d~ 
sirr0 J sm d 7 



«y 



Lastly, let A = */(l — sin 3 * sin s 7) , 

and in accordance with, the usual notation of elliptic integrals let 



F — 



y dy 



E = 



A $7. 



* « * * 



^0 



(5) 



Then we have the following transformations ;• 
cfo* 2 



*\ 



^ = 



qABG a sin 7 

du 2 



F 



sm^7 



a<ia> JqJABO "* a 8 sin 3 7 
d^r f 00 cfe ?2 



c?7 






" 
Sim 5 7 



a 8 sin 3 « sin 8 7 



(F-E) 



}- (6) 



5cZ& J JJi 8 C r ""a 3 sin 3 7j A 8 



% 



efo* 



2 



ccic ~"~J ABG & ~~* » 8 sin 8 7 



r tan 3 7 



€?7 



It remains to rednce the last two of (6) to elliptic integrals. 
If h and V be the modulus and its complement, the following are 
known transformations in the theory of elliptic functions, viz. .*— 



P*' flfy 1 -p & 3 sin 7 cos 7 

J fl S = P Pa ' 







tan 3 7 7 A tan 7 1 •« 
1 7 #* A/ 3 



e # e * 



HI olli >/ . J. 

ence | — -^ -$7=— 

A 8 7^ 



1-A S 







A 8 



dl ~~k*k'*® 



sm 7 cos 7 






(7) 



(8) 



(9) 
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In the present case &=sin a, &'=cos a, A=cos /3. Thus (8) and (9) 
enable us to complete the required transformation to elliptic in. 
tegrals of (6). 

Substituting then from (6) (8) (9) in the expression 

V=|m<* + — T- + -y- ir + — ,» >. 
L a da b do 



where 



c dc } 

m = f ?ra&c = |7ra 3 cos /3 cos 7, we have 

Tr « 2 ^ , 2 f # 2 ,_ _ N , / sin 7 cos 7 

asm 7 tt^in°7L sm * \cos~<xc< 

~j^v-V--^-(E--tan7C0S^)l . (10) 
;in 2 acosW cos 2 * ' J 



cos fi sin 2 * 



sin* 



Now suppose the ellipsoid to be rotating about the axis of z witli 
an angular velocity w, and let us choose the axes a, a cos <3, a cos 7, 
and the angular velocity w, so that the surface may be a surface of 
equilibrium. 

For this purpose V-J-|w 2 (a; 2 +^)=constant, must be identical with 



05 2 



t/2 sj2 

"o? + a 8 cos 8 /8 + a 2 cos 1 7 



= 1. 



Now in (10) we have V in the form 

V=Lx* + My*+Nz* + P, . . . 

whence a\L -f | w 2 ) = a 2 ( Jf + 1 &> 2 ) cos 2 0= a 2 JV cos V 

Hence Jv-ilf+JVcos 2 7tan 2 £=0, "| 

J^ 2 = N cos 2 7 — 1/ , 
or ^tt> 2 sin 2 /3=U cos 2 £— I/. 

There are two kinds of solutions of these equations (12). 
First, since 



(11) 



(12) 



d* 
da 



Jj—wbc-r- = — wa 3 cos i3 cos 7 . ,» . tT 

cos 8 cos 7 i v sin 2 7 ., 

= — 2?r £- — M -—d^ 

sin°7 Jo A 



Sltl'fy 

a°sm°7j0 A 



* 



3r=7rac-7r = — wa 3 cos 



^ ^ 2 f Y sin 2 7 



> • (13) 



cos fl cos 7 

— /Lit r— r x 

ein'V/ 



sm"7 



8in^7 7 
t-^7, 







A* 



it is obvious that Zr— ilf vanishes when *=0. 
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And since when x vanishes, /3 also vanishes, the equation 



L -M -f N cos 2 7 tan 2 /3=0 



is satisfied by 



a = 0, 



/3=0. 



That is to say there is a solution of the problem which makes 

Thus there is solution which gives us an ellipsoid of revolution. 
When sina=0, we have also /3=0, A = l, and 



o 



L= 



_ Z7r cos 7 



sin°7 



• 07 t cos 7 . . . 

sin- 7 dy — - -.--—. — ( sin 7 cos 7 — 7 ) , 

SlIl'V/ 



N 



2wcos7f 7 9TC0S7 /o ot \ 

= — ,- tan 3 7 dv = — — .— (27—2 tan 7) . 

sin ,5 7 J sm-^7 



Therefore 



|w 2 =^/"cos 2 7— L, 



7T 



= ta^^ 2<Y_2 tan 7 ~ ^ + tan2 7)( sm 7 cos 7-7)]» 



7T 

1 — 5- F°/ (3 -f tan 2 7) — 3 tan 7I, 
tan' 1 7 L " IJ7 



(14)* 



and the eccentricity of the ellipsoid of revolution is sin 7. 



To find the other solution when a, is not zero, we have by com- 
parison between (10) and (11), 



-\ 



L . - -7_ = E-F, 

27T COS /3 COS 7 

sin 2 a sin 3 7 sin 2 a sin 7 cos 7 ^, _ 
-^ • o 7, —= ~ — '—T-^ + F-Esec 2 ^ 

^7rCOS/^OOS7 COS 2 aCOS/3 

>T sin 2 a sin 8 7 n .„ 

iV . 77 =tan 2 a(E--tan7eos/3). 

2zr cos /3 cos 7 v if/ 

Hence the first of (12) gives 

(2F-E)+Esec 2 * 

1 i„ o 1 o,-, /-n i ^\ sin 2 a sin7 cos 7 

+tan^tan 2 /3cos 2 7(E— tan 7 cos /3)- ' ' 



> • 05) 



COS 2 <x COS ft 



0, 



* Compare with Thomson and Tait's ' Nat. Phil./ § 771, (3) ; or any other wori 
which gives a solution of the problem. 
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or 

E sec 3 *[l + (sin * tan p cos 7 ) 3 ] - (2E-E) 

~sec 3 *sin*tan^eos7(l + sin 3 /3)=0. . . . (16) 

In order to adapt this for computation, we may introduce the 
auxiliary angles denned by 

tan g==sin*tan/3cos7, tan S=am.p, . . . . (17) 

and the equation becomes 

E sec 3 * sec 3 £-- (2E-E) -sec 3 * sec 2 £ tan£=0. . . (18) 

The second of (12) gives 

£? . ^*™tl -- tan»*ooB 3 7 (B-taii7008i3)-(E-F), 
4?r cos p cos 7 i\ i rj \ j> 

, w 3 E—E Ecos 2 7 cos 8 cos 7 

whence -j = -r-5 — r-^- + — r-* 



w 3 



*7T 



47rcos^cos7 sin 2 asin 3 7 sin 3 7C0S 3 * cos 2 *sin 3 7 

=cot /3 cosec /3 cot 7(F— E) -j-cot 3 7 cos p sec 3 *E — cos 3 /3 cot 3 7 sec 2 *. 

. . . (19) 



Some of the subsequent computations were, however, actually 
made from a formula deduced from the third of (12), which leads to 



£i> 3 



A = cot /3 cot 7 cosec 3 /3(l -fcos 3 /?)(E— E) — cot 3 /3 cot 7 tan 3 * cosec /3E 

+ cot 3 /3 cot 3 7 sec 3 *. . . (20) 

By subtracting (20) from (19) we can deduce (16) ; hence it follows 
that (19) and (20) lead to identical results. Most of the subsequent 
results were computed from both (19) and (20), thus verifying the 
solution of (18). 

The formulae (18) and (19) are suitable for finding the solution, 
except when cc is small or nearly 90°, when the elliptic integrals 
become awkward to use. I have, therefore, found approximate 
formulae for these cases, but as the algebraic process necessary to 
establish them is somewhat tedious, the details are given in a note.* 



* Approximate Solutions of the Problem. 
From (7) we have 
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Approximate Solutions of the Problem — continued. 
Now, since (16) may be written 



J- -»-» •^-«-« -»~«N «* sin y cos y tc » _ 
pjE- (2F-B) j^—Jl^ — gm*y oos*y 



k*A 



A 2 



1 _ A , 

— E 4- • — ■ tan y 

tf2 ^^ Jc'2 r 



it follows from (a) and (8) that it may be written 



ry . 



A 2 



jl (5fy = sin^y COS 2 7 

A 3 



tan 2 y , 



(») 



o 



Again, we may write (19) thus :■ 



a; 2 Ac0S7 /T , _. -Acos'Vi-i- A 2 cos 2 7 

— _. ~(F— E) h -E- — ; - 

/tw ^ 2 sin 3 v ' & /2 sin 3 7 & /2 sin' J 7 



4ht ^ 2 sin 3 7 

A cos 7 
sin 3 7 

_ A cos 7 
sin 3 7 



-Y . 



sm 2 7 A cos 3 7 tan 2 y 



o 

r-rv 



A 



-<^7 



sin 2 7 
™A~ 



sin 3 7 J A 

n 



t^7 



<$7 — cos : 



„ [ tan s 



7 



(^7 



#©*>*& 



(«) 



Tt is easier to develop the equations when written in the forms (5) and (c) than 
when we work directly from the elliptic integrals. 

Write for brevity 

&='snra, g— cos <x, _p = cos7, £ = siny, Q=tany, A=log e cot(i7r — £7). 

First, when 7c is small— 

The following definite integrals are required : — • 



1 . »-l 

q n dy^ — — q n ~ i pH 

-- n n 



q n ~ 2 dy. . 



« e 



<<0 






n— 1 



C?7- 



(»-l) U^- S *y. 



a © # s © $--*.- 



(«) 



1 1 8 3 5 

Erom (d) and from — =^—j—=i + -tf^ + ^4tf + ■ • • 



A 2 



_ — coy sss . o ,ii « «— ■ — 

A 3 r AT F 4.2 



A 2 = l™-h 2 q 2 i and 
* 4 1.3 



o.JL o 

42 ? + 2 



fl'i' + Tsr + o* 2 



1 5 5 , 

6 * ^ 6.4 



5.3 
' 6.4.2^ + 6 



.3.1 \ 
— 71+. 

.4.2 / 



(/) 



£L^__^-^ + g 7+ *. 



S 6J2~4/ 



2 5 P 



k 6\4k2~4JS 



2 3 i> 



3.5.3 / 3.5.3. 1 3.1 

2.6.4,2 ^ + 7 ^2.6.4. 



2 ~4.2 r 



+ . . 



1 00 



-ie^-^p+y 




3 < 

•—2" 

8* 



! )J 



1 1 

Again, — — 1 + -&V+ . . »> and by (e) 



W 
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Approximate Solutions of the Problem — continued. 

Jt^-j^-J - y +\* [j !+ l« , -t T ] +• • • 

i O 2 1 r 3 3 ^ "1 

(1 3 3 \ 
~4;M 5 + 4 } P?-^yp 2 <l i )+ -. (*) 

The equation (5) is^V i^^=A 2 U» whence, equating (^) and (A), 
5 , 3 3 /, 8 , A 79 r , , 11 ,15 /15 3,3 Al 



5 3 /3 \ 

If A? be zero, we have -<f<p + -pg -7 1 -- +2?V ) = °- 

4 8 \o / 



+ . . .=0. 



m ,. ., , , sin 27— T ^sin4y ,. x 

This easily reduces to 7= — - . — — - — —t (*) 

J 1 — £ cos 4y ' 

of which the solution is 7=54° 21' 27", as stated in the text. 
Now 

ch\ J^ + jjM-rfg +i>YJ = s -2[ 1 - 2cos2 7 + cos ' l! r + 'ysin 4y], 

and with the above value of 7 is equal to —01355014. 
Also, with this value of 7 the coefficient of & is 0*0160432 ; so that 

00160432sin (y~54° 21' 27") =0'1355014& 2 , 

or, sin 2 «=10- 9266528 sin (y-5*° 21' 27 ,/ ), 

which is the equation (21) of the text. 



f 2 1 11^1 q q\ 

Again JX d7= __ 2p + - T+ -^_ 4 -2» p _L ?p+ _T,j, 

l/l 3 3 3 \ 



Therefore 



1^-^-^.(1-4^-4^-^.(^-1^]. 



Now As=l — JA: 2 2 2 — . .. 
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Approximate Solutions of the Problem — continued. 



Therefore 



^{-|q+(| + |q^]4^[>-tQ + (^ Q2+?2Q2 ) 7 ]- * w 

Hence from (c) and (y) 

which is the equation (22) of the text. 

Secondly, let Jc be nearly unity and g small — 

Then we require the following definite integrals : — 



[ 2^1 dy -JL aQ n-l *L [? n —dy 

"I 



Ldy = — q + A. 
V 



(Jc) 



(0 



Now 



1 1 fl 3 of \ 



<? 4 <7 4 ^ _ 9 7 



'd~ ~3 



5 -r,i+". 



By (*) and (0— 



p.i 2^ ]) 






Now 
therefore 



A 2 



i> 



.=1 +/Q 2 , 



A 2 f «? * T 3 [" 1^ 9 /15 13 ^„ 1 ^A"l 3 r « /IS V 



But 
therefore 



1*' 



-=1+Q 2 , 



|j> - 1,^^^ 



28 41 1 Y 



8 



24 



-| A [l + Q^r(f + fQ^Q^)] (») 
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Again, 



1_ 1 

A 



Approximate Solutions of the Problem— 'continued. 
±~ p* P~2 S \P* P 3 ) 

+ (> + |,-1a)]... 



The equation (b) for determining the axes is 

3 
Hence equating (wt) and (n) and dividing by r # we have 

which is the equation (23) in the text. 
The equation (e) for n> may be written 



(») 



= 0, 



«2 A 

4*™ Q 3 _ 



(i + Q 2 )|-f<*7 



fM?4 



Now 



A 



r a . i „ri . i i_.~ 



(i+Q 5 ) 



-dy 



-{(H 



a+QO + a** 



|a(i + q^)-|(i + qv]} 






VOL. XLI. 



Z 
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If a bo infinitely small, so that the Jacobian ellipsoid of three 
unequal axes becomes in the limit an ellipsoid of revolution, we have 
7 given by 

_ sin 2<y — j\ sin 4<y 
l—|eos4<y 



7 = 



The solution of this is 7 =54° 21' 27". 

If we write tan 7 = /, as in Thomson and Tait's ' Natural Phi- 
losophy,' §778', this equation becomes 

tan- 1 / ^ 1 + -W 2 

/ ""1+yfH-f 4 ' 

which is the equation (9), §778', of that work. 

The ellipsoid of revolution of which the eccentricity is sin 54° 21' 27" 
belongs to the revolutional series of figures of equilibrium, and is 
the starting point of the Jacobian series of figures. As shown by 
Sir William Thomson, it is the flattest revolutional figure which is 
dynamically stable. The Jacobian figures of equilibrium are initially 
stable, and as stated by M. Poincare,* there is for this value of 7 a 
crossing point of the two series, and an exchange of stabilities. 

If a. be small, it appears that sin a, is given by 



sin 2 * = 10*9266528 sin ( 7 _54° 21' 27"), 



(21) 



and w by 

,2 



OJ 



4t _ - i cot 3 7 [ (3 + tan 2 7 ) 7-8 tan 7] : 



bw 



sin 2 * 



+ ^t^y C( * sin27 ~" )tan7+ ' y( -» ? ~ ( »~ sin ^ )tan ^ )] • (22) 



Now 
Therefore 



(i + Q 2 ) 



Hence 



Approximate Solutions of the Problem — continued. 

a=*>(i+^ 2 qO. 



~dy 



-I 



Q 2 



dy 



A/3 



+ 2' 






which is the equation (24) in the text. 
* " Sur r^quilibre d'une Masse do Fluidc, &c," l Acta Mathematical 7, 1886. 
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These formulae are, it must be admitted, of but little use, since it 
would be necessary to take in higher powers of sin 3 # to obtain results 
for a variation of 7 of more than 1°. 

If a be near 90°, so that cos & is small, the approximate equation 
between <% and 7 is 

1 + tan 3 7 -j~cos 2 <x( \ 5 + y tan 3 7 +•£ tan 4 7 +-§• tan 6 7> 

: — log e cot(l7r— 17) . [H-|tan 3 7 + cos 3 a( J 8 5 -+-V-fcan 3 7 + tan 4 7)] 

sin 7 l 

=0. . . . (23) 

And to is given by 

u —L— — 11 __ — log^o^ljr—lry) . (|+tan 3 7)—f sec 2 7 



4nr 2 tan 3 7 
1 



-f--J-COS 3 <% 



SU17 



logtfCot (\rr — 17) . (f + 2 tan 3 7 + tan 4 7) — |— - y- tan 3 7 



i tan 4 r 



J . . (24) 



We shall return later to a modification of (24) which will be 
applicable to very long ellipsoids of equilibrium. 

Besides the angalar velocity and the axes of the ellipsoid, the other 
important functions are the moment of momentum, the kinetic energy 
of rotation, and the intrinsic energy of the mass. In order to express 
these numerically we must adopt a unit of length, and it will be con- 
venient to take a, where 

a 3 = abe = a B cos /3 cos 7. 

Thus a = a(sec /3 sec 7)*. 

Let g be the density of the fluid which has hitherto been treated as 
unity, and let (f 7r<r)£a 5 /i,, (-§• 7nr) 2 a 5 e be the moment of momentum and 
kinetic energy, then 

(%7r(r)%& 5 /j.=±m(a2 + W)iv=^ E 7rca?(sec p sec 7)^(1 + cos 2 /3)(47r<r) 4 ( JL_ ) , 

\4i7T(T/ " 

Thus /* = -|V3(sec/3 sec ^(1 + cosW^LY. . . . (25) 

The function (25) is the quantity which will be tabulated. 

,2\-l 



Again (f 7nr) 3 a 5 <?=^(f tt<t) V| /* . w = J v^dwff^a 5 . /* ( -^~) 

\4<7r(T/ 



> 



/ w 2 \ 4 

so that € =2V 3« I - — I . . . . . (26) 

z 2 
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Tlie function (26) is the quantity which will be tabulated. 

Thus in the tables the unit of moment of momentum is taken as 
(f7nr)*a 5 , or mta*, and the unit of energy as (|7ro-) 2 a 5 or ra 2 /a. 

It remains to evaluate the intrinsic energy, or the energy required 
to expand the ellipsoid against its own gravitation, into a condition 
of infinite dispersion. 

If di be an element of volume, then this energy is 



1 

2 



Yadt 



integrated throughout the ellipsoid. 

This will be denoted by (f7rer) 2 a 5 (i— 1), or m*sT l (i— 1), so that i 
will be positive. 

Now Y=Lx 2 -\-My 2 -\-Nz*-\-P, and if we denote by A, B, 0, the prin- 
cipal moments of inertia of the ellipsoid, we have 



1 



^W£=KB-fC-A)=!ma 3 , 



and similarly, 



Also 



y*(T dt~j;mb%, 



z*<rdt=±mc®. 



%) nJ •) 



(rdt~m. 



j j j 



m 



Hence — (i - 1) = T \m [La* + Mb* + Nc* + 5P] 



a 



= T Vma 2 (sec ft sec ^[L +M cos 2 £-f Ncos^-f- 5Pa~ 2 ], 
But if we take the values of X, M, N given in (15), and note that 



P = ttoPcos ft cos 7 . 



-F, 



a sin 7 
it easily follows that 

£ + Jf cos 2 /3 + iVeos 2 7+Pa- 3 =0. 
(i— -1) = — |wa 2 (sec/5sec7)^ . Pa" 2 

= — 1 ma"- (sec /3 sec 7)? . § -— 3 



Hence 



rn* 



a 



-F 



a^ a sm 7 



3 w- 2 (cos ft cos 7)*™ 

sm.7 



a 



Therefore 



% 



(cos /* cos 7)"-™ 
sin 7 



(27) 



For a sphere 7 becomes infinitely small, and F becomes equal to 7, 
so that F/sin7=l. Thus i — 1=— f. Therefore the exhaustion of 
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energy of a sphere of radius a is % m 3 /a ; which is the known result. 
For an ellipsoid of revolution a =0, and /3=0, and F = 7; so that 

The function (27) is the quantity tabulated below. It seemed pre- 
ferable to tabulate a positive quantity, and it is on this account that 
the intrinsic energy corresponding to the infinitely long ellipsoid is 
entered as unity. 

Having now obtained all the necessary formulas, we may proceed 
to consider the solution of the problem. 

We have to solve 

sec 3 a sec s £ , E-(2F-E)-tanfsec 2 «sec 2 a=:0, . . (28) 

where tan £= sin a tan /3 cos 7, tan B = sin ft = sin a sin 7, 

f 7 d P 

and F= ~, E= cosfldfy. 

JoCos/3 Jo 

The axes of the ellipsoid are 

- = (sec a sec 7)3, - = -cosi3, - = -0087. • (29) 

a v ^ Iy a a aa ' v/ 

If e 1? e 2 , e 3 are the eccentricities of the sections through c&, cb t ah 
respectively, we have 

e^sin/3, e 3 =sin7, e 3 = cos a sin 7 sec ^3. • (30) 
Having obtained the solution, we have to compute 

cot /3 cosec /3 cot <y(F— E) +cot 3 /3cos /3 sec 2 «B— cos 3 /3 cot 3 7 sec 3 «* 

. . . (31) 



A.7TG 



Then we next compute p and e and * from the formulas (25), (26), 
(27). 

The functions F and E are tabulated in Table IX of the second 
volume of Legendre's ' Traite des Fonetions Elliptiques/ in a table of 
double entry for a and 7 for each degree. 

The solution of (28) by trial and error was laborious, as it was 
necessary to work with all the accuracy attainable with logarithms of 
seven figures. 

The method adopted was to choose an arbitrary value of 7, and 

* As stated above, some of the computations were actually made from the 
formula (20). 
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then by trial and error to find two values of * one degree apart, one 
of which made the left-hand side of (28) positive, and the other 
negative. 

The smallest valne of 7 is 54° 21' 27", but after that value integral 
degrees for 7 were always chosen. 

The solutions for 7=55° and 57° could not be found very exactly 
froni the elliptic integrals with logarithms of only seven figures, but 
tlie solutions were confirmed by the approximate formula? (21) and 
(22). The solution for 7=80° was confirmed by the approximate 
formulae (23) and (24), and that for 7=85° was only computed 
therefrom, since when 7=80° the approximate formula gave nearly 
identical results with the exact one. 

The solution obtained is embodied in the table opposite. The first 
three columns give the auxiliary angles 7, «, f$, from which the remain- 
ing'results are computed. 

As a graphical result is much more intelligible than a numerical 
one, I have given two figures, showing the three principal sections in 
two cases, namely, where 7=60°, and 7=75°. For these figures a is 
taken as 2 cm., so that the volume of fluid is -f ?r x 2 3 cubic cm. 
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It will be noticed that tlie longer the ellipsoid the slower it rotates. 
It is interesting to observe that while the angular velocity continually 
diminishes, the moment of momentum continually increases. The 
long ellipsoids are very nearly ellipsoids of revolution about an axis 
perpendicular to that of rotation. Thus in fig. 2 the section through 
b and c is not much flattened. 

The most remarkable point is that there is a maximum of kinetic 
energy when a/a is about 3, or when the length of the ellipsoid is 
about five times its diameter. However, notwithstanding this maxi- 
mum of kinetic energy, the total energy always increases with the 
length of the ellipsoid. 

The kinetic energy is the product of two factors, one of which 
always increases, and the other of which always diminishes ; thus it is 
obvious that it must have a maximum. The result was, however, 
quite unforeseen, and it seems worth while to obtain simpler formulas 
for the case of the long ellipsoids. This may be done by taking as the 
parameter a/a, or the length of the ellipsoid, instead of ty. 

From the table we see that in the later entries /3 is very nearly 
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equal to 7, and that « becomes very nearly equal to 90°. Hence we 
may put «=90°, and /3=7. 
Thus, approximately, 

- = (sec j3 sec 7)*== (sec 7)* 
a 

and °«"»=(s)*' '<=*"-($■ 

The axes of the ellipsoid are 

Now if in formula (24) we only retain the higher powers of tan 7, 
we have 

a* 2 sin 7 cos 7 f It ,/, 1 \ *T 

— = — _L 1 __ log e cot (J*.-^) — J- 

47r<r tan 7 L.s in 7 J 



But 



cos <v r ~i 

= * "im 7 [* l0geC0t (^-27) -sin 7 J 



log e cot(|7r~i7)==log e -^^^==log e 2-f|log e ~. 

cos 7 a 



rn- 



Therefore writing 1— flog e 2=C, so that 0=0*3573, we have 

2 



4zr<7 



=*©'h-H 



If we put a/a= 5*042, this formula gives a> 2 /47nr= 0*01264. The full 
value in the preceding tables was 00131 ; thus even with so short an 
ellipsoid as this, the results agree within 4 per cent. With rougher 
approximation we have 



" ,3 o / a Y\ a 
- — =s-gl — 1 log,--, 



of which the limit, when a is large, is zero. 
For the moment of momentum we have 



3^ / «3 \4 

/* =y(sec p sec 7 )t(l + cos 2 /3) \^J 

=llr© ! 04)(><-«)' 
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or, with rougher approximation, 



3 AWi <*V 



of which the limit is infinite. 
Again, 



~~2 ^Kfrro) 



=*;['4]h-H- 

Now the function -( log e C ) has a maximum, when log e -=l + C 

ct \^ a / a 



1-3573, that is when -=1G96. 

a 



On comparison with our tables it is obvious that the approximation 
is bad, and that the true solution for a maximum is considerably dif- 
ferent from the above. Nevertheless this investigation shows that 
there is actually a maximum of kinetic energy. 



Since 
we have 



Y= log,cot Q*-h) =* [log|+f log,2] , 
If we like we may express these several results in terms of the 



a* 



minor and major axes of the ellipsoid, for b=c=- , and therefore 

a* 



a 3 =c%. 



Thus 
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